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文章 编号 :1005-3085(2010)04-0621-06 


图 的 最 大 二 等 分 问题 的 秩 二 松弛 算法 的 改进 


张 ” 芳 ， 徐 成 贤 
(西安 交通 大 学 理学 院 ， 西 安 710049) 


摘 要 : 本 文 在 吸取 半 定 规划 松弛 和 秩 二 松弛 方法 的 优点 ， 克 服 其 缺点 的 基础 上 ， 针 对 模型 目标 函数 非 凸 
的 特点 ， 提 出 了 图 的 最 大 二 等 分 问题 的 秩 二 松弛 模型 。 由 于 该 模型 变量 的 数目 没有 增加 ， 因 此 该 
方法 对 求解 大 规模 问题 很 有 优势 。 数 值 实验 表明 ， 这 种 算法 无 论 是 与 半 定 规划 松弛 还 是 原 秩 二 松 
弛 算法 相 比 ， 在 获得 目标 函数 值 相 当 的 情况 下 ， 运 行 时 间 较 短 。 

关键 词 : 图 的 最 大 二 等 分 问题 ， 秩 二 松弛 ; 拟 Newton 法 

分 类 号 : AMS(2000) 90C27; 90C59 中 图 分 类 号 : 0221.2 文献 标识 码 : A 


给 定 无 向 图 G = (V, E) MEREXI REE W = (wsij)nxn， 其 中 7 = {1,2,… ,n} 为 节点 集 
&, ENWRS, WE (ij) 瓦 ， 定 义 在 此 边 上 的 非 负 权 值 为 wy = wi BU wy = 0。 图 
的 最 大 二 等 分 问题 就 是 把 节点 分 成 个 数 相同 的 两 部 分 9 和 3 = V\S (这 里 nn 是 偶数 )， 使 得 两 
节点 分 别 在 3S，5 中 的 边 的 权 值 wj 的 总 和 最 大 。 引 入 二 值 变量 zi, i = 1,2,… ,n。 若 i € 5， 
Way; =1; Hi E S, Wri = 一 上 。 Wa = (£1, T2, aas ee R” 表示 元 素 全 为 1 的 列 向 量 ， 
则 图 的 最 大 二 等 分 问题 可 以 表述 为 


1 
w* = max 》 qual — 223) 
i,j 


n 
st. Soa =0, m € {1,1}, i=1,2,---,n. (1) 
i=l 


图 的 最 大 二 等 分 问题 属于 组 合 优化 问题 ， 它 在 组 合 优化 中 起 着 非常 重要 作用 ， 在 网 络 优 
化 ， 工 程 技术 中 也 有 着 广泛 的 应 用 背景 ， 许 多 网 络 问题 都 可 以 转化 为 图 的 最 大 二 等 分 问题 ， 
例如 超大 规模 集成 电路 (VLSI 的 设计 技术 等 。 由 于 该 问题 很 难 获得 精确 的 解 ， 是 NP-hard 问 
题 。 近 年 来 ， 人 们 已 经 不 满足 常规 算法 求解 复杂 问题 ， 提 出 了 各 种 启发 式 近 似 算法 ， 其 中 松 驰 
算法 就 是 非常 重要 的 一 种 。 建 立 在 半 定 规划 (SDP) 松弛 基础 上 的 随机 近似 算法 可 以 在 多 项 式 时 
间 内 获得 原 问 题 的 一 个 好 的 上 界 ， 该 算法 对 于 图 的 节点 不 太 多 时 表现 得 很 有 效 ， 但 是 随 着 图 的 
节点 数目 的 增加 ， 转 化 所 得 的 半 定 规划 松弛 模型 变量 的 数目 呈 几 何 级 数 增加 ， 这 给 问题 的 求解 
带 来 了 难度 。 

Burer 等 中 提出 了 最 大 制 问题 的 秩 二 松弛 算法 ， 其 中 求解 图 的 最 大 割 问 题 的 秩 二 松弛 模型 时 
使 用 的 是 最 速 下 降 法 ， 理 论 分 析 和 实际 计算 表明 ， 最 速 下 降 法 的 效果 是 很 不 理想 的 。 本 文 吸 取 
半 定 规划 松弛 及 秩 二 松弛 方法 的 优点 ， 克 服 其 缺点 的 基础 上 ， 提 出 图 的 最 大 二 等 分 问题 的 秩 二 
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松弛 方法 ， 由 于 图 的 最 大 二 等 分 问题 的 秩 二 松弛 模型 变量 的 数目 没有 增加 ， 所 以 这 种 松弛 方法 
对 求解 大 规模 问题 就 很 有 优势 。 但 是 该 模型 导致 目标 函数 非 凸 ， 本 文 针 对 该 模型 的 特性 ， 提 出 
一 种 改进 的 求解 方法 。 数 值 实验 表明 ， 这 种 方法 是 有 效 的 ， 与 半 定 规划 松弛 方法 相 比 ， 获 得 的 
原 问 题 的 近似 最 优 值 相 近 的 情况 下 ， 运 行 时 间 较 短 。 与 原 秩 二 松弛 算法 相 比 ， 改 进 的 秩 二 松弛 
算法 获得 目标 函数 值 相当 的 情况 下 ， 运 行 时 间 较 短 。 








2 图 的 最 大 二 等 分 问题 的 秩 二 松弛 模型 及 特性 


首先 给 出 图 的 最 大 二 等 分 问题 的 秩 二 松弛 模型 。 我 们 用 单位 向 量 w e RRS) 中 的 
标量 zi， 用 向 量 内 积 vz ui 代替 标量 zizj， 于 是 限制 条 件 € {一 1,1} 就 变 为 |lvijl。 = 1， 其 
中 juil 表示 向 量 w 的 2- 范 数 。 这 样 ， 所 有 的 向 量 凡 都 在 R? 的 单位 贺 上 。 由 于 wi 是 二 维 单位 
向 量 ， 故 可 以 在 极 坐标 系 下 设 


cos 6; S 
vi = ， i=1,2,-- n. 
sin 6; 


Fe IK n ARAE, ,vn E R, WUE Rh n HE = (01,02,… ,0n)， 而 
且 对 任意 的 i, j= 1,2,---,n, v2 vj = cos(b 一 0;)。 若 定义 对 称 和 矩阵 Y (0) : Y(0) = (cos(0; 一 
0;))nxn， 则 对 任意 的 9€ R? 可 以 定义 


f(0) = SPI cos(0; — 8;), (2) 
i=1 j=l 
MBIA (eTax)? = (eeT) o (zx7) = eT Xe 二 0 成 立 ， HP X = raf = (zizj)nxn， 用 矩阵 Y(9) 代 
OX, TERRA ED a WHA 


=, oh 
min z (9) 
s.t. eTY(0)e=0, (3) 


该 松弛 模型 可 以 看 做 是 用 秩 二 限制 rank(Y) < 2 代替 半 定 规划 松弛 模型 中 的 秩 一 限制 ， 因 此 称 
为 图 的 最 大 二 等 分 问题 的 秩 二 松弛 模型 。 该 模型 是 一 个 具有 nn 个 变量 的 约束 最 优化 问题 ， 也 就 
是 说 变量 的 数目 没有 增加 ， 这 正 是 该 松弛 方法 的 优越 之 处 。 

下 面 分 析 秩 二 松弛 模型 的 特性 。 函 数 f (0) 的 二 阶 偏 导数 为 


a2 (8) Wij cos(0; — 0;), 若 ; FJ; 
80,00; ) — >> we cos(Ok — 0j), isj. 
k#j 
TRA 
G(6) = V2f(0) = W oY (6) — Diag( [W o Y (#)]e), (4) 


Hp “o” RE Hadamard 内 积 ， 即 矩阵 相应 元 素 的 乘积 形成 的 矩阵 ，Diag(z) 表示 由 向 量 z 的 
分 量 作为 对 角 元 素 形 成 的 对 角 和 矩阵 。 由 (4) 式 可 知 目 标 函 数 f(9) 的 Hessian ER G(9) 奇异 ， 所 
以 图 的 最 大 二 等 分 问题 的 秩 二 松弛 模型 的 目标 函数 非 凸 ， 不 能 保证 获得 原 问 题 最 优 值 的 一 个 
界 ， 为 此 本 文通 过 增加 迭代 次 数 和 对 初始 值 的 扰动 来 提高 目标 函数 的 值 。 
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3 ” 秩 二 松弛 算法 


Goemans 和 Williamson 在 文献 [6] 中 求解 图 的 最 大 割 问题 时 提出 了 随机 化 方法 ， 而 且 他 们 
还 证 明 : 如 果 图 的 所 有 权 值 非 负 ， 用 此 随机 化 方法 产生 的 割 值 的 数学 期 望 至 少 是 最 大 割 值 
的 0.87856 倍 。 本 文 仍 采用 Goemans 和 Williamson 的 随机 化 方法 。 假 定 我 们 从 模型 (3) 中 获得 
局 部 (或 全 局 ) 极 小 点 09， 定义 相应 于 6 的 二 等 分 值 记 为 


WO) = > Twy (1 = cos(0; a 0;)). (5) 





由 三 角 函 数 的 周期 性 知 ， 对 于 9 的 每 一 个 分 量 一 致 旋转 一 个 角度 r， 函 数 f(9) 的 值 不 变 ， 
即 对 任意 的 + € R, f(0)=fO+ re)， 也 就 是 说 F(6) 是 关于 每 一 个 分 量 b 的 周期 函数 ， 周 期 
为 2r。 一 般 地 ， 可 设 bi e (0,27), i=1,2,--- ,n(n 是 偶数 )， 且 901 < 0s。<…… < 0,,。 我 们 取 整 
Hk € [1,n/2), > 


-1 1, #ie[k,k+n/2), é 


-1, 其它. 
这 样 得 到 的 z = (21,22 tn) 满足 图 的 二 等 分 条 件 ， 记 相应 的 二 等 分 的 值 为 


r(x) = D Twyll — Ti£j), (7) 
ij 
为 了 十 分 有 效 地 取 遍 (6) 式 所 有 可 能 值 x:， 并 找 出 最 好 的 二 等 分 ， 则 可 以 利用 下 述 程序 来 实 
现 


步骤 1 HAO c R”, 满足 0; € [0,27), i = 1,2,---,n, HAO < bg < < On 
令 MB=0,k=1; 

步骤 2 由 (6) 产 生 一 个 等 分 制 z， 并 由 (7) 计算 r(z); 

步骤 3 若 r(z) > MB, MWS MB =r(z), x* =x, BI), HIR 4: 

步骤 4 若 k > n/2 -1， 则 停止 计算 ， 否 则 天 = 大 十 1， 转 步骤 2。 

从 上 述 程序 可 以 看 出 来 ， 由 程序 产生 的 二 等 分 仅仅 依赖 于 BR 中 单位 贺 上 的 点 的 顺序 ， 而 不 
是 这 些 点 的 实际 位 置 。 模 型 (3) 中 的 约束 只 是 对 点 的 位 置 限制 ， 而 对 他 们 的 顺序 没有 起 到 约束 
作用 ， 因 此 没有 必要 求解 约束 松弛 模型 (3)， 而 直接 求解 去 掉 约 束 条 件 的 松弛 模型 ， 这 样 仍然 
有 可 能 产生 相同 或 更 好 的 二 等 分 ， 并 且 花 费 的 时 间 要 短 。 

按照 Goemans 和 Williamson 的 理论 ， 由 程序 得 到 的 二 等 分 值 r(z*) 至 少 是 秩 二 松弛 模型 的 
最 优 解 9 对 应 的 二 等 分 值 (9) 的 0.87856 倍 ， 即 7(z*) > 0.87856w(0)。 由 于 不 能 保证 (0) 是 最 
大 二 等 分 值 的 一 个 上 界 ， 无 法 得 到 算法 具体 的 近似 比 。 但 是 在 某 种 程度 上 可 以 保证 ， 我 们 获得 
的 (9) 的 局 部 极 大 (或 者 说 f(0) 的 局 部 极 小 ) 越 好 ， 产 生 的 二 等 分 值 就 可 能 越 好 。 

KH Newton 法 求解 无 约束 优化 问题 需要 计算 目标 函数 的 Hessian 矩阵 Gy, 和 求解 线性 方程 
组 V2zF(zto)d = 一 Vf(z(*))， 计 算 量 大 。 而 拟 Newton 法 中 BFGS 校 正 公 式 是 迄今 最 好 的 
拟 Newton 公式 ， 对 于 一 般 的 函数 ， 它 都 具有 很 好 的 性 质 。 它 是 由 f(z) 的 函数 值 和 一 阶 导 数 
值 构造 对 称 和 矩阵 By, 来 近似 代替 Hessian E Gr, HDI de) = —Bi lgk 来 确定 搜索 方向 ， 避 
免 了 Hessian 阵 的 计算 。 为 了 减少 工作 量 ， 避 免 求解 线性 方程 组 do = -Bg 来 确定 搜索 
Aid, WA, = Be, Wd = -Hrg WA: = Hy + OMe, AH, 是 修正 和 矩阵， 并 要 
R Ay ys WEH Newton 条件 Hry ® = 8s, FRBBRKF Hr 的 BFGS 公式 
YOT Hy OT Hy OT + (Hpy®s®T)T 

JOTA) ) rr TFT » 68) 





Ags, = Hk + (1+ 
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设 Ay, EWIE, yOTs > 0， 则 由 式 (8) 得 Anas 仍 是 正定 矩阵 ， 且 当初 始 矩 阵 Ho E 
定时 ， 对 一 切 的 k&，KHi 是 正定 和 矩阵， 从 而 保证 了 拟 Newton 方 向 d%) 总 是 下 降 方向 。 于 是 
有 BFGS 拟 Newton 法 : 

步骤 1 HAs, Ho =I, k=0; 

步骤 2 如果 满 足 某 终止 条 件 ， 则 输出 z( 中 ， 停 止 ， 否 则 do = 一 有 ok 

步骤 3 ”进行 线性 搜索 获得 步 长 ok， 令 ze+D = 2) + and, 

步骤 4 thy =VF(ct) —VF(e), s = ht) — gl), 

步骤 5 由 (8) 式 计算 Hiri; 

步骤 6 上 = 大 十 1， 转 步骤 2。 

大 量 的 数值 实验 表明 BFGS 方法 与 低 精 度 的 线性 搜索 方法 结合 使 用 ， 能 收 到 较 好 的 计算 结 
果 。 因 此 在 确定 了 搜索 方向 之 后 ， 本 文选 择 非 精确 线性 搜索 Armijo 试探 法 站 确定 步 长 ax。 这 
样 我 们 就 可 以 根据 秩 二 松弛 模型 (3)，BFGS 拟 Newton 法 及 上 述 程序 建立 改进 的 图 的 最 大 二 等 
分 问题 的 启发 式 近似 算法 : 

步骤 1 MAN, NEERA, 0 < R", 1 = 1 MB=0; 

步骤 2 Wo 为 初始 点 ， 用 BFGS 方 法 解 去 掉 约 束 条 件 的 松弛 模型 (3)， 得 到 f(9) 的 极 小 
点 0; 

步骤 3 ”由 程序 计算 一 个 关于 6 的 图 的 最 好 的 二 等 分 z， 由 (7) 计算 r(z); 

步骤 4 若 r(z) > MB, WS MB =r(z), =r, BM, BPRS; 

步骤 5 若 1 > N， 则 停止 ， 否则 随机 取 A0€ R", F0=04+A0, 1 =141, BHR3. 


4 ”数值 实验 
本 文采 用 BFGS 方法 求解 去 掉 约 束 条 件 的 松弛 模型 (3)， 得 到 S(O) 的 极 小 点 ， 取 终止 条 件 为 
IV F(A) < 10-4， |Ja*tY) — 6] < 10-4(1 + 06®))), 


取 Ho 为 单位 矩阵 。 采 用 非 精 确 线 性 搜索 Armijo 试探 法 确定 步 长 ， 取 6 = 0.75, p= 02, y= 
6。 将 算法 用 Matlab 6.5 编程 ， 取 9° 为 随机 向 量 。 做 如 下 数值 实验 ，w 表示 近似 最 优 解 ，Cpu(s) 
表示 CPU 运行 时 间 ( 秒 )， 是 迭代 次 数 。 

首先 考虑 文献 中 出 现 的 几 个 问题 ， 问 题 1: 文献 四 中 的 一 个 图 。 问 题 2: 文献 [3 中 的 一 
个 图 。 问 题 3: 图 为 20 个 点 19 条 边 的 一 条 路 径 。 问 题 4: 10 个 节点 的 两 个 完全 图 之 间 用 一 
条 边 连 接 形成 的 图 。 问 题 5: 一 个 六 角 星 图 (节点 1 与 节点 2,5,6 连 接 ; 节点 2 与 3,6 连 接 ; 节 
点 3 与 4,6 连 接 ; 节点 4 与 5,6 连 接 ; 节点 5 与 6 连接 )。 问 题 6: 文献 四 中 的 一 个 图 。 对 于 这 
几 个 问题 ， 如 果 节 点 间 有 边 连接 ， 则 权 值 为 1， 否 则 为 0。 由 于 版 面 的 限制 ， 本 文 只 给 出 参考 
文献 ， 不 再 另 附 表 说 明 。 数 值 实 验 结果 见 表 1。 从 表 1 中 可 以 看 出 ， 对 于 图 的 最 大 二 等 分 问 
题 ， 除 了 问题 3， 用 秩 二 松弛 算法 得 到 的 近似 最 优 值 都 要 好 ， 运 行 时 间 一 致 较 短 。 

其 次 选取 W 为 随机 生成 的 稀疏 权 值 对 称 矩 阵 ， 并 用 2 代替 非 零 元 素 ， 将 秩 二 松弛 算法 与 改 
进 的 算法 作 比较 ， 运 行 结果 如 表 2。 从 表 2 中 可 以 看 出 ， 秩 二 松弛 算法 与 改进 的 秩 二 松弛 算法 
相 比 ， 节 点 数 为 200 时 秩 二 松弛 算 用 了 相当 长 的 时 间 ， 运 行 过 程 中 出 现 条 件数 较 大 的 警告 ， 但 
改进 算法 运行 时 间 并 不 长 ; 总 体 上 说 ， 节 点 数 为 400 时 得 到 的 目标 函数 值 稍 微 差 一 点 ， 但 是 运 
行 时 间 一 致 较 短 ， 而 且 对 于 图 的 节点 数 越 多 ， 改 进 算法 的 效果 越 明显 。 
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表 1: 文献 中 的 几 个 问题 的 数值 实验 结果 





秩 二 松弛 算法 (N = 1) 改进 的 秩 二 松弛 算法 (N = 8) 























表 2: 秩 二 松弛 算法 与 改进 的 秩 二 松弛 算法 的 比较 








图 秩 二 松弛 算法 (N = 1) 改进 的 秩 二 松弛 算法 (N = 1) 
w Cpu(s) 
0. 160 


n w Cpu(s) 









































































10 46 0.731 

20 172 1.342 0.080 
30 402 2.984 0.040 
40 704 2.835 0.061 
50 16.804 0.180 
60 18.257 0.181 
80 82.418 0.330 
100 144.178 0.441 
200 1110.239 1.873 
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由 于 图 的 最 大 二 等 分 问题 的 秩 二 松弛 模型 变量 的 数目 没有 增加 ， 所 以 这 种 松弛 方法 对 求解 
大 规模 问题 很 有 优势 ， 但 是 秩 二 松弛 模型 导致 目标 函数 非 凸 ， 不 能 保证 获得 原 问题 最 优 值 的 一 
个 上 界 ， 为 此 通过 增加 迭代 次 数 和 对 初始 值 的 扰动 来 改善 目标 函数 的 值 。 数 值 实验 结果 表明 ， 
对 相同 的 测试 问题 ， 与 半 定 规划 松弛 算法 相 比 ， 用 秩 二 松弛 算法 获得 的 目标 函数 值 相 差 不 大 ， 
但 是 该 算法 运行 时 间 比 半 定 松弛 算法 要 短 ; 与 原 秩 二 松弛 算法 相 比 ， 改 进 的 算法 获得 目标 函数 
值 相当 的 情况 下 ， 运 行 时 间 较 短 。 
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Improvement on the Rank-two Relaxation Algorithm for the Graph 
Max-bisection Problem 


ZHANG Fang, XU Cheng-xian 


(School of Science, Xi’an Jiaotong University, Xi'an 710049) 


Abstract: A rank-two relaxation model is proposed for the graph max-bisection problem. The ap- 
proach is based on theories of SDP and rank-two relaxation approaches, and the approach utilizes the 
character that the objective function is non-convex in the model. Since the number of variables does 
not increase, this approach is suitable for large-scale problems. Extensive numerical experiments show 
that, compared with the SDP and the original rank-two relaxation approaches, our algorithm need less 
time to approach the similar value of objective function. 
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